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We consider the possibility of a large phase of Bs-Bs mixing in supersymmetric SU(5) and SO(IO) 
models. We find that in the SU(5) model, the magnitude of this phase is correlated with the 
branching ratio of r — > fi-y and the phase can be within 2a of the recent UTfit analysis. In the 
case of SO(IO) models, this correlation can be relaxed and a large phase can be obtained. In this 
scenario, a non-zero value of CP asymmetry for B — > Xs'y will be predicted. We predict the sparticle 
mass ranges for the LHC for these models once the UTfit result is accommodated and discuss the 
dark matter and the anomalous magnetic moment constraints on this analysis. 

PACS numbers: 12.10.Dm, 12.60.Jv, 12.15.Ff 



Recently, CDF and D0 collaborations have announced 
the analysis of the flavor-tagged Bs J /'^4> decay. The 
decay width difference and the mixing induced CP violat- 
ing phase, 0s, were extracted from their analysis 1]. In 
the Standard Model (SM), the CP violating phase is pre- 
dicted to be small, -/.^ = 2(3, = 2 arg (-14^^/^/^,^^*^ ~ 
0.04. However, the measurements of the phase are large: 



0s (D0) = 



[0.32,2.82] (68% CL) (1) 
0.57t°-jO(stat)t;5:EJ?(syst) (2) 



The UTfit group made a combined data analysis includ- 
ing the semileptonic asymmetry in the Bs decay, and find 
that the CP violating phase deviates more than ia from 
the SM prediction l2l| . This implies the existence of new 
physics (NP) and that the NP model requires a flavor 
violation in h-s transition. 

Supersymmetry (SUSY) is the most attractive candi- 
date to construct the NP models. In SUSY models, the 
flavor universality is often assumed in squark and slep- 
ton mass matrices to avoid large flavor changing neutral 
currents (FCNC) in the meson mixings and the lepton 
flavor violations (LEV) S] . Even if we assume the flavor 
universality, the non-universality is generated from the 
evolution of renormalization group equations (RGE) . In 
the minimal extension of SUSY standard model (MSSM) 
with right-handed neutrinos, the induced ECNCs from 
RGE effects are not large in the quark sector, while siz- 
able effects can be generated in the lepton sector due to 
the large neutrino mixing angles [3] ■ In the grand unified 
theories (GUT), the loop effects due to the large neu- 
trino mixings can induce sizable effects in the quark sec- 
tor since GUT scale particles can propagate in the loops. 
The patterns of the induced ECNCs highly depend on 
the unification scenario, and therefore, it is important to 
investigate the ECNC effects to obtain a footprint of the 
GUT models. 

If the quark-lepton unification is manifested in GUT 
models, the flavor violation in b-s transition can be re- 
sponsible for the large atmospheric neutrino mixing 5|, 
and thus, it has to be related to the r — > decay [6|. 
The branching ratio of the t f^j is being measured at 



the i?-factory, and thus, the future results of LEV and 
the phase of Bg-Bg mixing will provide an important in- 
formation to probe the GUT scale physics. In this Letter, 
we study the correlation between Br(T 11^) and (j), in 
SU(5) and SO(IG) GUT models, and investigate the con- 
straints in these models from the observations in order 
to decipher GUT models. 

We first describe the SU(5) and SO(IO) GUT mod- 
els which we investigate in this Letter. In the SU(5) 
model, the right-handed down-type quarks {D'^) and left- 
handed lepton doublets (L) are unified in 5 represen- 
tation. The quark doublets (Q), right-handed up- type 
quarks {W), and right-handed charged-leptons {E'^) are 
unified in 10, and the right-handed neutrino {N'^) is 
a singlet under SU(5). The superpotential which in- 
volves the Yukawa interaction is: Wy — Y^^ lOilO jH^ + 
yj-'lOiSj-ffg + Yj^^5iNjH5, where i,j denote the genera- 
tion indices, and H5 and are the Higgs fields in which 
colored Higgs fields are unified with the Higgs doublets. 
The charged-lepton Yukawa coupling Yg is unified to the 
down- type quark Yukawa coupling Yd, Ye — Yj, in the 
minimal SU(5) setup. In the basis where Ye and Majo- 
rana right-handed neutrino mass matrix Mjv are diago- 
nal, the Dirac neutrino Yukawa coupling is denoted as 
Y„ = UlY^^^suT^ = 1^ the unitary matrix 

Ul is same as a mixing matrix of neutrino oscillation, 
and it contains large mixing angles. The Dirac neu- 
trino Yukawa interaction can generate off-diagonal ele- 
ments of the (squared) scalar mass matrix for 5, M|, 
due to a loop diagram in which N'^ and propa- 
gate. The induced fiavor violating term in the scalar 
mass matrix is proportional to Yi,Yj , and M| can be 

parameterized as M? = m|[l — KC/Ldiag(fci, ^2, 
where 7715 is a scalar mass for 5 at a cutoff scale, 
M*, which is possibly the Planck scale, and k repre- 
sents the size of flavor violation in the scalar mass and 
K ~ (Y^'''^)l^/{8TT^){3 + Al/ml) InM./McuT in the uni- 
versal SUSY breaking scenario where mo is a universal 
scalar mass and Aq is a universal trilinear coupling. As- 
suming the Dirac neutrino Yukawa coupling is hierar- 
chical (fci,fc2 ^ 1), we obtain the 23 element of M| 



2 



as — 1 /2 TTig K sin 2023 e*" , where 6*23 is a 2-3 mixing in 
Ul, which is large and responsible for the large atmo- 
spheric neutrino mixing unless there exists a fine-tuned 
relation among Y^^^^ and Mj^. The phase e*" generates a 
phase of SUSY contribution for Bg-Bs mixing amplitude, 
Mf^^^, and the absolute value ofMf^^^ is controlled by 
K sin 2023- The off-diagonal elements of the scalar mass 
matrix for 10 are also generated by colored Higgs loop, in 
addition to the MSSM contribution, but they are small 
since they get generated from quark mixings. 

In the SO(IO) GUT model, which we use in this Let- 
ter, all matter species are unified in the spinor repre- 
sentation 16 Q. The matter representation can couple 
to the 10, 126 and 120 Higgs representations, Wy = 
h.jWaQjW + f^jlQil6jT26 + h[jW,Wjl20. The Ma- 
jorana neutrino mass is generated from the / coupling 
term when the B — L direction of 126 gets a vacuum ex- 
pectation value (VEV). The fermion Yukawa couplings, 
Yu,d,e,u, are given by the linear combinations of h, /, 
h' multiplied by Higgs mixings. If there is no cancella- 
tion among the h, /, h' couplings, the Dirac neutrino 
Yukawa coupling does not have large mixing in the ba- 
sis where the charge-lepton mass matrix is diagonal, due 
to the presence of right-handed neutrino in 16. Thus, 
the large neutrino mixings may not originate from YJ^. 
However, they can originate from the relative mixing an- 
gle of h and / couplings in the model. Such a structure 
can be constructed even in the SU(5) models with addi- 
tional Hig^ fields. The SO (10) model, however, is more 
predictive [8| since the Majorana coupling and the con- 
tributions to Dirac fermion masses are unified to the / 
coupling. In the basis where Ye is diagonal, the / cou- 
pling is denoted as U f'^^^^U'^ , and the unitary matrix 
U is the neutrino mixing matrix when the triplet term 
of the type II seesaw [9| is dominant. Through the loop 
diagram in which the heavy particles from 126 (or 120) 
propagate, flavor violations term can be generated for all 
sfermion mass matrices and it is proportional to //^ (or 
h'h''^). The off-diagonal elements are similarly denoted in 
the SU(5) case. The phase of SUSY contribution for the 
Bg-Bs mixing amplitude Mf2^^^ can be generated from 
the phases of 23 elements of A/? and M^^. The two 
phases are independent in the SO(IO) model in the basis 
where Yd is real diagonal, and the freedom of choosing 
one of these phases governs the phase of Mfg^^^. It is 
important that the absolute value of MfJ-"^^ is enhanced 
if both left- and right-handed squark mass matrices have 
off-diagonal elements [1, [13] ■ Therefore, a large phase of 
Bs-Bs mixing can be expected in SO(IO) model which is 
much larger compared to the SU(5) model. 

Now we discuss the phase of Bg-Bg mixing in the GUT 
models. We use the model-independent parameteriza- 
tion of the NP contribution: Cs.e^^'^^- = Mlf^/Mf^, 
where 'full' means the SM plus NP contribution, Mfj^' = 



ized by two real parameters Cb^ and (fiB^- The time 
dependent CP asymmetry {S = sin 0s) in Bg — s- J/ipcj) 
is dictated by the argument of A/^j" '■ <t>s — — argMf2'\ 
and thus (ps, is the deviation from the SM prediction: 
(f)s — 2{Ps — (pBs)- It is important to note that the 
large SUSY contribution is still allowed even though the 
mass difference of Bg-Bg ll| is fairly consistent with the 



SM prediction. This is because the mass difference can 
be just twice the absolute value of M^"^. The consis- 
tency of the mass difference just means Cb^ ^ 1, and 
a large 4>Bs is still allowed. Actually, when Cb^ — 1, 
the phase is related as 2sin(^B^ ~ A^^ /J^^, where 



A 



NP.SM 



NP.SMi 



In the model-independent global 



Mf/. The NP contribution can be parameter- 



analysis by the UTflt group, the fit result is 

AY/Af^ e [0.24, 1.38] U [1.50, 2.47] (3) 

at 95% probability [2]. The argument of being free 
in GUT models is due to the phase in off-diagonal ele- 
ments in SUSY breaking mass matrix (in the basis where 
Yd is a real diagonal matrix), and one can choose an ap- 
propriate value for the new phase in the NP contribution. 
Therefore, the experimental data constrain A^^ /A^^, 
and that means k sin 2^23 is constrained for a given SUSY 
particle spectrum. 

Among the ffavor violating decay modes, the current 
bound for Br(r — > /X7) is less than 4.5 x 10~^ [13] and 
Br(6 ^ 57) = (3.55 ± 0.26) x 10"'' In fact, in the 

models where gaugino masses are unified at the GUT 
scale (neglecting anomaly mediated SUSY breaking con- 
tribution), the T — > /i7 constraint is stronger than the 
b ^ constraint in most of the parameter space. This 
is because the squark masses are raised much more com- 
pared to the slepton masses due to the gaugino loop con- 
tribution since the gluino is heavier than the Bino and 
the Wino at low energy. This gaugino effect is also im- 
portant to allow a large phase in the Bg-Bg mixing. The 
gaugino loop effects are ffavor invisible and they enhance 
the diagonal elements of the scalar mass matrices while 
keeping the off-diagonal elements unchanged. If the fla- 
vor universal scalar masses at the cutoff scale (7715 , mio in 
our notation) become larger, both Br(r — > ^7) and (pB, 
are suppressed. However, Br(r ^j) is much more sup- 
pressed compared to (pB, for a given k sin 26*23 because 
the low energy slepton masses are sensitive to TO5 and 
mio while squark masses are not so sensitive due to the 
gluino loop contribution to their masses. 

In the SU(5) model, we flnd that ms has to be larger 
than 1.2 TeV in order to obtain a large (j)B^ at 95% prob- 
ability of the UTflt result if mio,fj, < 1 TeV (/i is the 
Higgsino mass) for a universal gaugino mass 7tt.i/2 — 300 
GeV and tan/3 = 10, which is a ratio of Higgs fields' 
VEVs. On the other hand, in the SO(IO) model, we find 
that mie (= = toiq) needs to be larger than 500 
GeV if < 1 TeV and mi/2 = 300 GeV and tan/3 = 10. 
The reason of a smaller scalar mass being allowed in the 
SO(IO) model is due to a left-right enhancement effect 
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in the box diagram. The constraint on masses based on 
UTfit resuh provides an important guidance to search for 
the SUSY particles at the LHC. Actually, in general, if 
quark-lepton unification is manifested in the GUT mod- 
els, the slepton masses needs to be heavy (especially in 
the SU(5) model) in order to suppress r ^ ^7 and to 
obtain a large phase of Bg-Bs mixing. 

The diagrams for r — > which can provide impor- 
tant effects are the chargino loop diagrams. This con- 
tribution can be suppressed if 11 and/or ms are large. If 
mio is small (which means that right-handed sleptons are 
light), the Bino diagram can contribute. In SU(5) model, 
the off-diagonal elements of the right-handed sleptons are 
small, but the Bino diagram can generate LFV through 
left-right mixing of sleptons. The amplitude of r — > iij 
is proportional to tan/3, while (j)B^ does not depend on 
tan/3 much. Therefore, in order to obtain a larger (/)_Bs, 
heavier SUSY particles and smaller tan/3 are favored. In 
Fig.l, we show the contour plot for A^^ /A^^ when the 
Br(T — > fij) bound is saturated. To draw the figure, we 
choose mi/2 — 300 GeV, ms = 2 TeV, and tan/3 = 10. 
We assume that SUSY breaking Higgs masses, m//„ and 
ruHa: are free to make to be a free parameter, but we 
assume to mjj„ = mjy^ at the GUT scale just for sim- 
plicity. In the figure, the blue shaded region shows la 
allowed region and the yellow shaded region is enclosed 
for 2(7 allowed region. The right-bottom area in the fig- 
ure is excluded since the stability condition of the Higgs 
potential is not satisfied. In the area, closer to the ex- 
cluded region, the mass of the charged Higgs boson is 
light and therefore, 6 ^ 57 becomes large. If fj, is large, 
the gluino contribution becomes large for right-handed 
operators of 6 — > 57 (which is often called Cm, Csr) due 
to the large left-right mixing for sbottom, and in this 
region, 6 — > 57 is more important rather than r — > /17. 
One can see that large A^^/A^^ > 0.38 (which is at 
68% probability) is allowed for a large value of ^. This 
is because the chargino contribution for r — > /Lt7 is sup- 
pressed in those area. The contours are curved down for 
small TOio because the Bino diagram starts contribut- 
ing. It is worth noting that the stau coannihilation re- 
gion (where the lightest stau and the lightest neutralino 
mass difference is 5-15 GeV [ll|) mio 100 GeV to 
satisfy the dark matter content is not allowed at 95% 
probability in this figure. To revive the stau coannihila- 
tion region, tan /3 ~ 5 is needed and then scalar trilinear 
coupling has to be chosen appropriately to satisfy the 
lightest Higgs mass bound. The Higgsino dark matter 
is not very favored (though it is allowed for tan /3 = 10 
in the 95% probability) since a small value of fi does 
not suppress the chargino contribution for t ^ fi'j. The 
stop coannihilation mechanism or yl-funnel region would 
be more suitable to explain the dark matter content if 
the quark-lepton unification is manifested and the ex- 
perimental constraint from fiavor violation is considered. 
A detailed discussion on dark matter constraint will be 
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FIG. 1: Contour plot for ^/Af^ when Br(r ^17) bound 
is saturated. The detail is given in the text. 

presented elsewhere. 

Unfortunately, in order to reduce the chargino contri- 
bution to r — !■ /i7, the sleptons must be heavy enough, 
and as a result we find that one cannot provide a so- 
lution for a muon g — 2 anomaly [iBj to obtain large 
<j)Bs if quark-lepton unification is manifested. Therefore, 
we should consider the possible breaking of the manifest 
quark-lepton unification. So, let us consider the possi- 
bility that K sin 2^23 is different between the squark and 
slepton sectors. To begin with, the 23 mixing angle 6123 
can be different between quarks and leptons because Yd 
and Ye may not be simultaneously diagonalized. In min- 
imal SU(5) model, and Y^ are unified, but it gives a 
wrong prediction for the 1st and 2nd generation masses, 
and we need a correction from 45 Higgs field or non- 
renormalizable interaction. Actually, there is a freedom 
to choose 6'23^'°" <^ 6*23 ^'"^ which is needed to relax the 
constraint arising from r iij. However, the motiva- 
tion to explain the large neutrino mixing is lost. This 
situation is same as in the SO (10) models. Thus, let us 
consider the possibility that k is different between quark 
and lepton sectors. Actually, it should be different since 
the flavor violation terms are generated from the loop dia- 
gram in which heavy particles propagate, and the heavy 
particles should split when GUT symmetry is broken. 
In SU(5) model, however, it always gives wrong direc- 
tion, i.e., kI'''"'^ < k'<=p'°". The RGE effect to generate 
the flavor violation survives till the right-handed Majo- 
rana mass scale for left-handed slepton, but it ends at 
the colored Higgs mass scale for the right-handed down- 
type squarks. To satisfy the nucleon decay bounds, the 
colored Higgs need to be heavier rather than the right- 
handed Majorana mass scale when is less than 0(1). 
Therefore, in a model, where the flavor violation origi- 
nates from the Dirac neutrino Yukawa coupling, a large 
(j)B^ is disfavored. In this sense, the Fig.l is drawn con- 
servatively. When the Majorana mass is 10"'^'' GeV, the 
colored Higgs mass is 10^^ GeV, and the cutoff scale is 
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10^^ GeV, one then obtains 2k1"'^'''^ ~ ^lepton^ ^^lat 
time, the contour values for A^^ /A^^ should be reduced 
roughly by half to saturate the t ^ bound and then 
II has to be larger than 1 TeV. 

In the SO(IO) model, 6'23's for quarks and leptons are 
not very different if there is no huge cancellation in the 
Yukawa couplings, /, and h' . However, k can be dif- 
ferent among the sfermion species, and it depends on 
the SO(IO) breaking vacua and the Higgs spectrum from 
126. If the mass of one of the Higgs representation is light 
compared to the SO(IO) breaking scale, the off-diagonal 
elements for some fermion species are generated. For 
example, when the SU(2)/j symmetry remains unbroken 
and (1,1, ±1) fields from 126, which are S\]{2)r Hig- 
gsinos, remain light, the off-diagonal elements for only 
right-handed slepton are generated. Obviously, these do 
not generate a large (j)B,- When (8, 2, ±1/2) fields of 
126 are light, only off-diagonal elements for squarks are 
generated while the same elements for sleptons are small. 
This is the proper way to generate large 4>b, without suf- 
fering from the r — > fij constraint. It is interesting to 
note that the light (8, 2, ±1/2) fields are good candidate 
to suppress the nucleon decay [l3| in these models. In 
this way, it is possible that the experimental measure- 
ments of large (pB^ and Br(T — > fij) for a given SUSY 
particle spectrum can constrain the GUT scale particle 
spectrum. Thus, more experimental data is very impor- 
tant to probe the GUT scale physics. 

If a suitable SO(IO) breaking vacuum is selected, the 
r — !■ fij and relation can be completely broken. If 
both left- and right-handed squarks have off-diagonal el- 
ements in this vacuum, the large (j)B^ can be easily ob- 
tained. In Fig. 2, we demonstrate the allowed region 
for different ktoq. We assume a universal scalar mass 
TTiQ = ms = mio = "mn,, = rnn^- The figure is drawn 
in the case of tan/3 — 10, though tan/3 is not an im- 
portant parameter in the plot. When mo is small, the 
trilinear scalar coupling Aq has to be large in order to 
make k large. If t ^ /i7 is suppressed by the choice of a 
vacuum, the b ^ constraint becomes more important 
especially for large tan /3. However, since the phase of the 
chargino contribution is free in the SO (10) boundary con- 
dition (actually it is independent of the phase of Affj^^^), 
experimentally allowed solutions can be found as long as 
the gluino contribution for Ctr, Csr is not very large. 
When the chargino contribution is large, the direct CP 
asymmetry for B — > Xgj [13] may become large. This is 
similar to the case where (j)B^ is large while Cb^ "^1- In 
Fig. 3, we plot the direct CP asymmetry and A^^/A™ 
for allowed Br(i3 Xg'y)- Due to the freedom of the off- 
diagonal elements of up-type squark mass matrices, the 
b ^ s"f constraint can be satisfied even in the usually ex- 
cluded region in the case of minimal supergravity model. 
Actually the parameter mi/2 = 200 GeV we choose to 
draw the Fig. 3 is not allowed in the minimal super- 
gravity for tan/3 = 40. Thus in this region, one needs 
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FIG. 3: Direct CP asymmetry oi B —* Xs7 decay vs 
AT /Af^. We choose = 200 GeV, mo = 250 GeV, 

and tan j3 = 40. 

off-diagonal elements of squark mass matrix, and thus, 
a non-zero value of A^^ is predicted. If Br(i3 — > Xs^) 
is fixed to a particular value, this plot is just a circle. 
The dark matter and the anomalous magnetic moment 
constraints are easily satisfied in this scenario. 

In this Letter, we have emphasized the importance of 
T /17 and (f)B^ correlation in GUT models, since they 
can be correlated directly by 23 mixing. The constraint 
from /i — > 67 may be also important, but this Br calcu- 
lation highly depends on the details of flavor structure 
which can have a freedom of cancellation. Therefore, we 
have not talked about the /i ^ 67 constraint. We refer to 
the Ref. [l3| for an analysis of flavor violation including 
the first generation. 

In conclusion, we have investigated the large phase 
of Bs-Bs mixing by comparing SU(5) and SO(IO) GUT 
models. The existence of (fiB^ in GUT models can tell us 
whether the flavor violation originates from Dirac neu- 
trino Yukawa coupling or Majorana coupling such as the 
126 Higgs coupling in the SO(IO) model. At present, 
the SO(IO) models are more preferred rather than the 
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Dirac neutrino induced flavor violation in SU(5) models. 
It is important to note that we can distinguish these two 
scenarios once more experimental data on Bg-Bg mixing 
phase and Br(T — > /i7) decay are available along with the 
data from the LHC. 

This work is supported in part by the DOE grant DE- 
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